UNCLASSIFIED 

..  413000 

DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DTIC  CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


I 


NOTICE:  When  government  or  other  dravings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  In  connection  vlth  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  Incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  In  any  vay 
supplied  the  said  dravlngB,  specifications,  or  other 
data  Is  not  to  be  regarded  by  implication  or  other- 
vise  as  In  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  Invention  that  may  In  any  vay  be  related 
thereto. 


63-  y-  ^ 


NORMAL  MODE  VIBRATIONS  OF  SYSTEMS  OF  . 
ELASTICALLY  CONNECTED  PARALLEL  BEAMS  ' 

by 

J.  M.  Seelig 
and 

W.  H.  Hoppmann  II 


Department  of  the  Army  -  Ordnance  Corps 
Contract  No.  DA-30-069-AMC-195(R) 


Department  of  Mechanics 
Rensselaer  Polytechnic  Institute 
Troy,  New  York 


July  1963 


swat#  Wft&Ff&F  S***  -  -  .  -.-  «...  -’.  rr 


\ 


NORMAL  MODE  VIBRATIONS  OF  SYSTEMS  OF 
ELASTICALLY  CONNECTED  PARALLEL  BEAMS 

by 

J.  M.  Seelig  and  W.  H.  Hoppmann  II 
Rensselaer  Polytechnic  Institute 
Troy,  New  York 

ABSTRACT 


The  report  presents  the  development  and  solution  of  the 
differential  equations  of  motion  of  a  system  of  n-elastically 
connected  parallel  beams.  Using  the  results  developed  for  the 
general  n-beam  system,  the  particular  case  of  a  two  beam  system 
was  analyzed  in  detail.  The  frequencies  and  associated  mode 
shapes  for  various  support  conditions  are  presented  explicitly. 


Vibration  experiments  were  performed  in  order  to  ascer¬ 
tain  the  degree  of  applicability  of  the  theory.  It  was  found 
that  at  least  up  to  the  eighth  mode,  reasonably  good  agreement 
obtains. 
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INTRODUCTION 

Flexural  vibration  of  the  single  bar  or  beam  has  been  ex¬ 
tensively  studied  by  many  Investigators.  The  validity  of  the 
Bernoulli- Euler  theory  has  been  studied  and  interesting  gener¬ 
alizations  of  that  useful  theory  have  been  developed.  Little 
work  has  been  done  on  structures  built  up  from  beams.  A  par¬ 
ticular  case  of  Interest  consists  of  a  parallel  system  of  thin 
beams  elastically  connected.  The  dynamics  of  such  structural 
units  is  of  Interest  particularly  because  of  their  possible  use 
in  various  areas  of  technology,  including  that  concerned  with 
space  platforms  for  use  in  orbital  problems. 


Preliminary  to  studies  of  dynamical  response  of  such 
systems  It  Is  necessary  to  determine  the  normal  mode  shapes 
and  their  associated  frequencies.  It  Is  the  purpose  of  thlB 
report  to  provide  explicitly  these  data  and  the  experimentally 
determined  basis  for  their  validity. 


SYSTEM  OP  DIFFERENTIAL  EQUATIONS  OF  VIBRATION  OF  n-BEAMS 
Using  Bemoulli-Euler  beam  theory,  the  system  of  differ¬ 
ential  equations  for  n-elastically  connected  beams  shown  in 


Figure  1  can  be  written: 

52w, 

Eli  - »r  +  Pi  -^5-  -  -  k^Wj^  -  w0) 


dx 


at* 


aV 


a2w. 


aV 


aV 


EI  - jr  +  p_  - -  -  k_  ,  (w  -  w  , ) 

n  a jT  n  at2^  ^-1  n  n-1 


(1) 


“i  *  ci  7^  -  -  ‘'i-iK  -  wi-i)  -  ■'iK  -  wi+i> 


(3) 


According  to  the  general  theory  of  elasticity,  the  free 

vibrations  of  such  a  system  are  harmonic,  so  we  may  write  for 

each  w„  : 
n 

wn  -  Xn{x)ei0it  (4) 

where  to  is  the  natural  frequency  and  X^x)  is  a  function 
of  the  space  variable  x  only. 

Using  Eq.  (4)  the  system  of  partial  differential  equa¬ 
tions  (l)-(3)  reduce  to  the  following  system  of  ordinary 
differential  equations 


Figure  1. •  System  of  n-Elastically 
Connected  Slender  Beams 


-  k1(X1  -  Xg) 


(5) 


pA  - 


•  t 


Pi®2*! 


ki-l^Xi  "  Xi-1^"  ki^Xi  "  Xi+l)  ^ 


d4X. 


El 


n 


n 


dx 


Si-l^  "  Xn-1^ 


(7) 


4nth 

n 

I 

n-1 


Solving  these  equations  Simultaneously,  we  obtain  a  single 
order  ordinary  differential  equation  on  say  X-^ 


d4nX, 

“n-JTIn 


U0X1 


(8) 


where  the  un  are  in  general  simple  functions  of  the  physical 
parameters  of  the  problem  and  the  frequency  d>,  Eq.  (8)  is 
obviously  an  ordinary  differential  with  constant  coefficients 
whose  solution  may  be  written  as  usual 

X-l  -  ce*x  .  (9) 

Substituting  Eq.  (9)  into  Eq.  (8)  gives  the  character¬ 
istic  equation 
n 

X  ^  ^  “  0  •  (10) 

n-0 
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! 
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I 

} 


If  we  let 

P  «  *4  (11) 


Eq.  (10)  becomes 

I  S  e"  -  0 

n-0 


(12) 


Solving  Eq.  (12)  for  the  0's  and  using  Eq.  (ll),  all 
the  Vs  are  determined.  The  solution  may  now  be  written 


Ivv  • 


(13) 


m-1 


In  the  solution  of  Eq.  (12),  if  any  of  the  0n's  are 
zero  or  are  repeated,  the  solution  (13)  to  the  differential 
equation  will  be  modified  acoording  to  the  well  known  methods 
for  ordinary  differential  equations  with  constant  coefficients 
for  those  particular  values  of  0n  . 


If  solution  (13)  la  then  substituted  into  the  system  of 
ordinary  differential  equations,  we  can  get  Xg  , . . . ,  X^^  as 
follows 


J  m  1 #2, . . . ,n 
1  *  1,2,.*., 4n 


(14) 


The  summation  convention  is  used  for  the  repeated  index  where 
the  c^j  are  linear  functions  of  Vs,  o’s  and  physical  r 

constants, 
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Any  elastic  rotational  or  elastic  support  boundary  condi¬ 
tions  may  be  investigated.  For  each  of  the  beams  in  the  n- 
beam  system,  there  will  be  four  associated  boundary  conditions. 

Applying  the  4n  homogeneous  boundary  conditions  we  get 
a  system  of  4n  homogeneous  simultaneous  equations  in  the 
constants  c^ 


lis  ‘‘i 


0 


1  m  1,..., 4n 
s  -  1, . , . ,4n 


(15) 


where  the  a1(J  are  linear  functions  of  X'b,  <o's  and  physi¬ 
cal  constants. 


In  order  to  have  non-trivial  solutions,  the  determinant 
of  the  coefficients  of  the  Cj.  must  vanish. 


I*i.l  ■  °  •  <16> 

The  expansion  of  determinant  (l6)  gives  the  frequency 
equation.  For  each  natural  frequency  the  c ^  can  be  deter¬ 
mined  in  terms  of  one  of  the  other  cA  from  Eq.  (15) •  The 
complete  solution  is  therefore 


eiot 


1  ®  1  p  •  •  #  j  4n 

j  «  1,  •  e  e,  n 


(17) 


-6- 


FREE  VIBRATIONS  OF  ELASTICAIX*Y 
CONNECTED  DOUBLE  BEAM  SYSTEMS 

As  a  preliminary  to  an  experimental  verification  of  the 

theory  and  to  illustrate  the  general  theory  of  n-beams,  the 

case  of  a  double  beam  system  consisting  of  two  identical  beams 

elastically  connected  will  be  studied  in  detail*  The  equations 

of  motion  then  are: 

d^w,  a2w, 

El - *  +  p - ±  -  -  k(w,  -  wp)  (18) 

ax4  at2  1  2 


a^w2  a2w 

ec — r  +  p — r  “  -  k(w9  -  w,)  (19) 

ax4  at2  2  1 


The  solutions  may  be  taken  in  the  following  form: 

.let  (20) 


wx  -  Xx(x)e3 


w2  -  Xg(x)e 


lot 


(21) 


Substituting  Eqs.  (20)  and  (21)  into  Eqs.  (18)  and  (19)  give 

(22) 


d4X, 


El  — r  -  “  -  k(xi  -  X?) 

dx 

dv 

B 1  -  P®2Xg  -  -  k(Xg  -  Xx) 


(23) 


Solving  Eq.  (22)  for  Xg  we  get 
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(24) 


S  mi  - 

k  dx*  k  1 


and  now  substituting  Xg  into  Eq.  (23)  we  get 


(K)2  +  (IS  .  SSI  2?£)  ^  +  <«£  .  S&jx,  -  o  (25) 


dx 


Letting 


El 


a  and 


CD 


b  ,  Eq.  (25)  becomes 


d8xl  ,  {2a  -  2ab\  ^l  .  b2  -  2b  Q 


(26) 


All  Ends  Simply  Supported 

Using  the  procedure  outlined  for  the  n-beam  case,  we  will 
now  treat  the  case  of  two  beams  being  simply  supported.  The 
boundary  conditions  for  this  case  are 


X^o)  - 

0 

x2(a) 

xj(0)-  - 

0 

XJU) 

X^o)  - 

0 

Xgd) 

x£(o)  - 

0 

x|(i) 

(27) 


Prom  the  boundary  conditions  (27)  and  using  expression  (24) 
we  finally  get  the  determinants!  equation 
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Solving  the  determlnantal  equation  (28)  the  frequencies 
and  associated  mode  shapes  are  shown  to  he: 

(a) 

o?  -  S  (IE)4  (29) 

P  * 


*2 


A' sin 


mrx 


A' sin 


rnrx 

i 


(30) 


(b) 

«2  -  S  (iE)4  +  2k  (31) 

p  i  p 


X,  -  A' sin  £25 

1  I 

(32) 

X~  »  -A'sin 

^  I 


Following  a  procedure  exactly  analogous  to  the  one  out* 
lined,  the  mode  shapes  and  their  associated  frequencies  for 
the  following  cases  have  been  determined. 


(a) 


All  Ends  Free 
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where 

cos  k^l  cosh  k^l  -  1 


A' 


£- 


osh  kj^x  +  cos  k^x 


osh  1^1  -  cos  1 ^1  sinh  1^1  -  sin  k^A 


sinh  k^x  +  sin  k^x^ 


*2  "  X1 


(b> 

a>2  .  ^4  El  +  2k 
P  P 

where 

cos  1^1  cosh  k^l  «  1 


c- 


osh  kjjX  +  cos  l^x  sinh  k^x  +  sin  k^x^ 
osh  k^l  -  cos  1^1  sinh  J^A  -  sin  k^l^ 


*2  “  "  X1 

As  another  possible  solution  each  bar  of  the  system  may 
vibrate  as  a  rigid  body  both  angularly  and  transversely  with 
frequency 


There  also  exists  the  trivial  solution  where  the  system 
translates  and  rotates  as  a  rigid  body  without  vibration, 


-11. 


All  Ends  Fixed 


(a) 


a,2  -  kj*  S 


where 


cos  1 cosh  kni  -  1 

rsinh  k^x  -  |)  sin  k^x  -  |) 


X1  -  *'  { - JTI 

sinh 


sln 


0>) 


0,2  .  .4  EI  +  2k 
P  P 


where 

cos  1^1  cosh  k^f  -  1 

..  . .  r8lnh  *!><*  -  4> 

X1  -  A’  1 — .  .  M 


sinh 


sin  k^(x  -  -|) 

I3T 


} 


-  x. 


-12- 


Double  Cantilever,  Both  Beams  Fixed  at  0 


(a) 


“2  *  kn1‘  f 


where 


cos  kni  cosh 


-  1 


rcosh  l^x  - 
1  \:osh  k_l  + 


osh  kyjX  -  cos  k^x  slnh  l^x  -  sin  kjjX'i 


1^1  +  cos  slnh  +  sin  1^1 


I 


(b) 


+  ilE 

P  P 


where 


cos  cosh  1^1 


-  1 


fCOSh  k^x  - 

1  Usosh  k„i  + 


osh  k^x  -  cos  kj^x  slnh  k^x  -  sin  k^X'j 


kj^i  +  cos  1^1  slnh  1^1  +  sin  1^1 


- Bin  VI 

+  sin  k  V 


-  -  X, 


-13- 


-14- 


NOW  MC  INOICATID  M  A  PROPORTION  or  LCNOTM  L  MEASURED  FROM  LEFT  END 
CLASTIC  CONNECTORS  BETWECN  SAKS  ARE  NOT  SHOWN  EN0 

WEWLAR  NATURAL  FREQUENCY  ,  OJ-CtfJ  *  ^  ,  RA^»«> 

Figure  2.  Modal  Shapes  and  c  •  youngs  modulus  ,  l */wf 

Natural  Frequencies  for  1  *  ****  MO****T  <*■  •*‘**r*  or  mam  cum  section,  «n* 

Some  Systems  of  Elastically  U*  LENGrTH  BCAMt ,N 

_  ,  €  *  *****  ****  WN«T  LENGTH  OF  WAN!  LS  WK'IlH* 

Connected  Double  Beams  u.  ‘ 

K-  KAODOLU&  OF  CLA&TVC  CONNECTORS ,  |_*/|N* 

*  COCFFlClfHTS  PROM  TAM  ASQvE  TM( 
•osscRirr  n  refers  to  move  --twiV 
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Due  to  this,  the  8x8  determinant  was  reduced  to  the  evalua¬ 
tion  of  two  4x4  determinants.  If  however,  this  symmetry  Is 
not  present,  the  zeros  of  the  8x8  determinant  must  be  found. 
As  an  example  the  following  case  Is  shown. 

The  boundary  conditions  for  one  bar  simply  supported, 
the  other  free,  are 

x1(o)  -  0  X1(i)  -  0 

X£(0)  -  0  X"(l)  -  0 

(33) 

x£(0)  -  0  x£U)  -  0 

X£'(0)  -  0  x»'(l)  -  0 

Using  relationship  (24)  and  the  boundary  conditions  the 
following  determlnantal  equation  Is  obtained. 
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The  zeros  of  determinant  (34)  are  now  found  giving  the 
natural  frequencies.  From  these  terms  the  mode  shape  asso¬ 
ciated  with  each  natural  frequency  can  be  determined.  Por  a 
specific  case  they  are  shown  In  Figure  3* 

VIBRATION  EXPERIMENTS  WITH  DOUBLE  BEAM  SYSTEM 

The  experimental  model  used  (See  Figure  4)  In  the  normal 
mode  vibration  experiments  was  supported  on  a  massive  rein¬ 
forced  concrete  base  3'  x  3'  x  10'.  The  beam  supports  were 
made  of  3/4"  steel  plates,  in  which  provision  was  made  for 
the  placement  of  pin  and  roller  supports.  The  pin  support 
consisted  of  an  actual  pin  placed  through  the  bar  at  its  mid 
height,  while  a  roller  support  was  provided  by  use  of  four 
roller  bearings  placed  on  the  lateral  boundaries  of  the  bar, 
so  as  to  provide  a  momentless  end  condition  which  would,  how¬ 
ever,  allow  axial  motion  of  the  bar.  Tracks  were  placed  above 
and  below  the  model  and  were  firmly  attached  to  the  steel 
support  plates  giving  added  rigidity  to  the  supporting 
structure . 

The  model  tested  consisted  of  two  1/2"  x  1/2"  x  39.9" 
parallel  bars,  elastically  connected  with  20  coil  springs 
placed  on  2  inch  centers.  The  bars  were  cut  from  standard 
cold  rolled  sections  and  had  a  Young's  modulus  close  to  29  x 
10^  psl  as  determined  by  test.  The  springs  were  2.0"  in 
length  and  had  a  spring  constant  of  40  pounds  per  inch.  They 
were  affixed  to  the  bars  using  Duco  cement. 
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Figure  4.  Experimental  Apparatus  for  Determination 
•  of  Frequencies  of  Vibration  and  Nodal  Patterns 


The  system  was  driven  by  an  electromagnet  supplied1  from 
an  oscillator.  The  output  of  the  oscillator  was  displayed  on 
the  upper  beam  of  a  Type  502  Tektronix  Dual  Beam  Oscilloscope. 

The  response  of  the  model  was  obtained  using  capacitive  type 

2 

pickups  placed  on  the  aforementioned  tracks.  This  signal 
was  displayed  on  the  lower  .beam  of  the  oscilloscope.  Resonant 
frequencies  were  determined  by  examination  of  the  capacitive 
pickup  response  as  the  frequency  spectrum  was  traversed.  A 
marked  increase  in  response  indicated  a  resonant  frequency. 

Due  to  the  pull-pull  characteristic  of  the  electromagnet,  the 
frequency  response  of  the  beam  system  was  twice  the  frequency 
set  on  the  oscillator. 

The  upper  beam  vertical  amplifier  was  then  switched  to 
the  horizontal  deflection  plates  and  a  figure  8  Lissajous 
figure  resulted  due  to  the  2:1  frequency  ratios.  The  capaci¬ 
tive  pickup  was  then  moved  on  ’-•he  track  until  the  figure  8 
horseshoe  degenerated  to  a  s^aight  line.  Continued  movement 
of  pickup  inverted  the  figure  8  horseshoe.  The  nodal  point 
was  located  when  the  straight  line  configuration  was  observed. 
Both  bars  were  checked  at  the  various  resonant  frequencies 
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and  all  nodal  points  were  recorded.  In  order  to  determine 
the  phasing  of  the  vibrating  bars,  the  responses  of  two  ca¬ 
pacitive  pickups,  placed  at  corresponding  locations  of  the 
bars,  were  displayed  on  the  upper  and  lower  oscilloscope 
beams.  A  difference  switch  on  the  Type  502  algebraically 
subtracted  the  lower  beam  input  from  the  upper  beam  Input. 

As  one  pickup  was  above  the  beam  system  and  the  other  below 
it,  the  response  when  two  points  of  the  system  were  In  phase 
resulted  in  a  zero  amplitude  (straight  line)  sine  wave,  while 
when  they  were  out  of  phase,  a  sine  wave  whose  amplitude  was 
the  sum  of  the  amplitudes  of  the  original  sine  wave  responses 
was  observed. 

The  relative  amplitudes  of  displacement  were  determined 
using  micrometer  barrels  fitted  with  conical  caps.  The  mi¬ 
crometers  were  electrically  connected  to  a  unit^  which,  when 
contact  was  made  by  a  micrometer  with  a  bar,  a  clicking 
switch  as  well  as  a  neon  light  would  be  activated.  Initial 
readings  on  the  micrometers  which  were  set  at  various  loca¬ 
tions  on  the  bars  were  made  with  the  system  at  rest.  The 
system  was  then  driven  at  a  resonant  frequency.  The  microme¬ 
ter  heads  were  then  made  to  just  contact  the  vibrating  bars 
and  additional  readings  were  made.  By  differencing  these 


^W.  H.  Hoppmann  II,  "Impulsive  Loads  on  Beams",  Proc.  Soc. 
Exp.  Stress  Analysis,  10,  No,  1,  p.  157,  (1932). 
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readings,  the  relative  amplitudes  of  displacement  were  deter¬ 
mined. 

The  mode  shapes  were  drawn  using  the  results  of  the  nodal 
point,  phasing*  and  relative  amplitude  measurements. 

Three  cases,  each  with  different  end  supports,  were 
tested.  These  were:  Case  I  -  all  four  boundaries  simply 
supported;  Case  II  -  all  boundaries  free;  Case  III  -  upper 
beam  simply  supported,  lower  beam  free.  For  Case  II,  light 
string  was  suspended  from  the  upper  track  and  was  looped  on 
the  upper  beam,  at  locations  corresponding  to  the  locations 
of  nodal  points  for  the  particular  resonant  frequency  being 
investigated,  as  support  for  the  system.  The  location  of  the 
string  aupports  at  the  nodal  points  did  not  change  the  free- 
free  boundary  conditions  and,  therefore,  did  not  introduce 
any  error. 

As  a  check  on  the  pin  and  roller  supports,  a  single  beam, 
simply  supported,  was  driven  and  its  resonant  frequencies 
were  determined.  These  values  were  checked  against  calculated 
theoretical  values.  The  values  were  within  a  few  percent  of 
one  another,  indicating  a  reasonable  approximation  to  simple 
support  conditions. 

NUMERICAL  EXAMPLES  AND  EXPERIMENTAL  RESULTS 

In  order  to  provide  a  basis  for  comparison  between  calcu*- 
lated  results  and  experimental  results,  the  double  beam  system 
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chosen  corresponded  to  those  used  as  models  In  the  experimental 
work. 


For  purposes  of  Illustration)  three  cases  were  Investi¬ 
gated.  These  were:  Case  X  -  all  four  boundaries  simply  sup¬ 
ported;  Case  I"  -  all  boundaries  free;  Case  III  -  upper  beam 
simply  supported,  lower  beam  free. 

Using  the  formulae  shown  in  Figure  2  and  Figure  3*  the 
frequencies  in  Tables  1,  2,  and  3  were  computed.  The  experi¬ 
mental  frequencies  shown  in  Tables  1,  2,  and  3  were  determined 
by  the  procedure  outlined  in  previous  section. 

Case  I  -  All  Four  Boundaries 
Simply  Supported 

Calc  Exp 

v.p.s, 

28  29 

79  79 

113  112 

135  137 

253  244 

264  261 

451  456 

460  490 

Table  1 
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Case  II  -  All  Boundaries  Free 


Calc 

Exp 

v.p, 

.8, 

65 

64 

98 

97 

178 

173 

192 

191 

350 

334 

358 

354 

578 

574 

583 

590 

Table  2 


Case  III  w  Upper  Beam  Simply 
Supported,  Lower  Beam  Free 


Calc 

Exp 

v.p. 

s. 

23 

23 

49 

48 

72 

70 

84 

84 

128 

126 

205 

182 

Table  3 


For  Case  III,  If  Figure  3  Is  used  for  combinations  of  I, 
EE  and  k,  the  same  values  of  c  will  be  found  as  for  the 
particular  case  investigated. 

The  location  of  nod*i  points  for  all  cases  Investigated 
checked  to  within  2%  of  their  calculated  values. 


CONCLUSIONS  AND  DISCUSSION 

Theory  has  been  developed  for  the  free  vibrations  of 
n-elastically  connected  parallel  beams.  As  a  consequence 
the  normal  mode  shapes  and  associated  natural  frequencies  can 
be  calculated  for  any  combination  of  homogeneous  boundary  con¬ 
ditions.  For  a  very  small  number  of  beams  the  calculations 
can  be  conducted  with  simple  computational  devices.  However, 
for  a  larger  number  of  beams  it  is  clear  that  high  speed  com¬ 
puters  will  be  required. 

For  the  double  beam  system  the  frequencies  and  normal 
mode  shapes  have  been  determined  for  several  interesting 
boundary  conditions.  Results  have  been  tabulated  and  compared 
with  those  obtained  from  experiments  on  steel  beams  coupled 
with  steel  springs. 

It  is  concluded  that  for  the  lower  modes,  at  least  up  to 
the  eighth,  the  agreement  between  theory  and  experiment  is 
very  good, 
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